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Abstract 

We study the statistics of wind-driven open ocean currents. Using the Ekman layer model for the inte- 
grated currents, we investigate, analytically and numerically, the relation between the wind distribution and 
its temporal correlations and the statistics of the open ocean cuiTcnts. We find that temporally long-range 
con^elated wind results in currents whose statistics is proportional to the wind-stress statistics. On the other 
hand, short-range coiTclated wind leads to Gaussian distributions of the current components, regardless of 
the stationary distribution of the winds, and therefore, to a Rayleigh distribution of the current amplitude 
if the wind stress is isotropic. An interesting result is the existence of an optimum in the amplitude of the 
ocean currents as a function of the correlation time of the wind stress. The results were validated using an 
oceanic general circulation model. 

PACS numbers: 92.60.Cc, 91.10.Vr, 88.05.Np, 92.60.Gn 
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Ocean currents are generated by local and remote forces/factors, including winds, tides, buoy- 
ancy fluxes and various types of waves. While many studies have investigated the distribution of 
the wind over land llHSD, focusing on its relevance to energy production, the distribution of ocean 
currents has received much less attention Moreover, currently there is no accepted theory 
explaining the observed statistics of surface ocean currents. Here, we propose a simple physical 
theory for the distribution of wind-driven ocean currents and its relation to the spatially variable 
temporal correlations of the wind (see Fig. [T]). We show that the distribution of wind-driven ocean 
currents strongly depends on the temporal correlations of the wind-when the wind exhibits long- 
range temporal correlations, the ocean current statistics is proportional to the wind-stress statistics 
while, for short-range correlations of the wind, the different components of the current vector fol- 
low Gaussian distributions. It was previously reported that the probability density function (PDF) 
of ocean currents follows the WeibuU distribution [|4j4Sl] — we argue that this is not necessarily the 
case, even if the wind stress magnitude follows the WeibuU distribution [|9|]. 
Oceans play a major and important role in the climate system, and ocean circulation underlies 
many climate phenomena from scales of meters to thousands of kilometers and from scales of 
minutes to decades. The gap in our understanding of ocean current statistics leads to a lack in the 
understanding of the statistics of other related climate variables. Filling in this gap will be useful 
in many fields: it may help to predict extreme current events and, thus, may help to securely de- 
sign maritime-associated structures. In light of the increasing efforts to find alternative sources of 
energy and the idea of using ocean currents as such a source, it is necessary to know the currents' 
PDF to better estimate the energy production and to appropriately design ocean current turbines 
that will withstand even extreme current events. Moreover, knowledge regarding current statis- 
tics and, in particular, its relation to the wind stress statistics may improve the parametrization of 
small-scale processes in state-of-the-art general circulation models (GCMs). 
The study of wind-driven ocean currents goes back more than 100 years, to the time when Ekman 
I lo| proposed his classical simple model to explain the effect of Earth's rotation on upper ocean 
currents. His model predicted that the depth integrated current vector is perpendicular to the wind 
vector, a prediction that was largely proven by observation. Since then, many studies have used 
Ekman's model to propose more realistic models for ocean currents as well as for surface winds 



12l1 . In what follows, we use Ekman's model [1^ to study the statistics of wind-driven ocean 



currents. 

The Model: In spite of the simplicity of the Ekman model, it will enable us to start investigating 



2 



Z3 



60°N 



40°N - 



20°N 



0° 



20°S - 



40°S 



60°S 




0° 100°E 160°W 

Lonqitude 



60°W 



I 



n 



40 

30 

20 

1 

9 

8 



FIG. 1 : A map of the correlation time (in days) of the wind stress magnitude. We define the correlation time 
as the time at which the normalized auto-con^elation function first drops below 1/e. Clearly, there is a large 

and 



vaiiability in the correlation time. One can also notice the remarkably shorter correlation time over 
compared with that over the ocean. The map is based on the NCEP/NCAR reanalysis six-hourly data il3l1 . 
Other definitions of the correlation time yield qualitatively the same results. The contour line indicates a 
one day correlation time. 

the coupling between the wind stress and the ocean currents. The time, t, and depth, z, depen- 
dent equations of the Ekman model lioll . describing the dynamics of the zonal (east- west) U and 
meridional (south-north) V components of the current vector, are: 



dt 
dV 

'dt 



(1) 



where / = {A-k/Tj) sin(0) is the local Coriolis parameter (T^ is the duration of a day in seconds 
and is the latitude), and u is the eddy viscosity coefficient, assumed to be homogeneous. The 
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boundary conditions are chosen such that the current derivative, with respect to the depth coordi- 
nate z, is proportional to the integrated current at the bottom of the layer described by our model 



and to the wind stress vector {t^, r„) at the surface IH 
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Here, we introduce the following notation: h is the depth of the upper ocean layer, r is a pro- 
portionality constant, (r^,, Ty) are the wind stress components, and po is the ocean water density 
(assumed to be constant hereafter (po = 1028A;5f/m'^)). By integrating eqs. ([T]) over a sufficiently 



deep layer (i.e., h ^ ^Ju/2f), we obtain the equations describing the integrated currents and their 
coupling to the wind stress. 



du 
dt 



fv — ru + 



Po 



— = -fu - rw H . 

dt Po 

To allow a simpler treatment of these equations, we define w 
obeys the following equation: 



dw 
'dt 



—ifw — rw + 



T 



Po 



(4) 



u + iv. It is easy to show that w 



(5) 



where r = r^. + iTy. The formal solution of equation ([5]) is 



1 

Po 



(6) 



This formal solution demonstrates that the currents depend on the history of the wind stress, and 
therefore, the distribution of the currents depends, not only on the wind stress distribution, but on 
all multi-time moments of the wind stress. However, the two extreme limits are quite intuitive. 
When the correlation time of the wind stress, T, is long (i.e., T ^ 1/r and T ^ 1//), one 
expects that the currents will be proportional to the wind stress since the ocean has enough time 
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to adjust to the wind and to almost reach a steady state. In terms of eq. in this limit, r(t') 
can be approximated by r(t) since its correlation time is longer than the period over which the 
exponential kernel is non-zero. In the other limit, when the correlation time of the wind stress is 
very short (i.e., rT <^ 1 and fT <^ 1), namely, the wind stress is frequently changing in a random 
way, the ocean is not able to gain any current magnitude. In this case, the central limit theorem 
implies that each component of the current vector is Gaussian distributed. 

To allow analytical treatment, we proceed by considering two special idealized cases in which 
the formal solution takes a closed analytical form. For simplicity, we only consider the case of 
statistically isotropic wind stress, namely, (r") = (r^') where is a positive integer and the 
different components of the wind stress are independent-(r™r^') = where m, n, are positive 
integers. 

Step-like wind stress: A simple way to model the temporal correlations of the wind is by 
assuming that the wind vector randomly changes every time period T while remaining constant 
between the "jumps". Integrating the solution of w (eq. Q) and taking into account the fact 
that the distribution of w at the initial time and final time should be identical, one can obtain the 
expressions for the second and fourth moments (ensemble average over many realizations of the 
stochastic wind stress) of the current amplitude. To better relate our results to the outcome of 
data analysis, we also need to consider the time average since the records are independent of the 
constant wind stress period. Calculating the double averaged second moment we get 

2\ I l~Ai r,r(l-Ai)+/A2 



iP + r^) Po 

where we have introduced the notations Ai = exp(— rT) cos(/T) and A2 = exp(— rT) sin(/T). 
The two extreme limits discussed earlier can be easily realized and understood. When the correla- 
tion time T is very long compared with 1/r, the second moment of the current is proportional to 
the second moment of the wind stress, namely, = (|t|^)/[(/^ + r'^)pl]. This is exactly what 

one would expect. The other limit is when rT, /T ^ 1 and, in this case, {\wf) ~ {\rf)T/{2rpl), 
which is very small and approaching zero as the correlation time approaches zero. This limit is 
also what one would expect for a rapidly changing wind that cannot drive significant currents. The 
double averaged fourth moment is 
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where D = exp(— 2rT) . While at the limit T — )• both the second and the fourth moments vanish, 
the ratio between them remains finite and is equal to 2. This corresponds to the Rayleigh distri- 
bution which originates at the fact that the components of the currents are independent and each 
has a Gaussian distribution with a zero mean and the same variance. One can easily understand 
the origin of the Gaussian distribution by considering the central limit theorem which can be ap- 
plied in this limit. Each of the current components is the integral of many independent, identically 
distributed random variables (the wind stress at different times). Based on the second and fourth 
moments in the limit of T — )■ 0, we conjecture that the overall PDF of the ocean current speed, in 
this case, is given by the Rayleigh distribution. 

Orenstein-Uhlenbeck wind stress: The second private case we consider is the Orenstein- 
Uhlenbeck wind stress. This process results in a Gaussian wind stress with an exponentially 
decaying temporal correlation function. It is important to note that, due to the Gaussian nature of 
the force, once we have the two-point correlation function and the first two moments of the wind 
stress, we have all the information on the driving force. The correlation function of the wind stress 
components is given by 

{n{t)r,{t')) = 5.,{t!) exp(-7,|t - t'\). (9) 
The long time limit of the current components' averages are: 

The long time limit of the current amplitude variance, = — (m)^ — )^, is given by 

I^yPlr [pi {r + i?) ■ ^^^^ 



Under the assumption of isotropic wind stress, the variance may be expressed as ~ 
(|rp) (r + 7) / (ppr (/^ + (r + 7)^))- Due to the Gaussian nature of the wind stress and the 
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FIG. 2: (a) The second moment of the cuiTent amphtude versus the wind stress con^elation time, T. Ana- 
lytical results (solid lines) compared with the numerical solution of the Ekman model (empty symbols) and 
the MITgcm modeling (blue symbols filled with either black or red) of the cun^ents in a simple artificial 
lake. The wind stress was drawn from a Weibull distribution with two different values of the k parameter 
as specified in the figure, (b) The Weibull k parameter of the current distribution versus the constant wind 
stress duration, (c) The Weibull k parameter of the current distribution versus the k parameter of the wind 
stress distribution. The Coriolis parameter, /, in this case, is 0, corresponding to its value at the equator, (d) 
Same as panel (c) but /, in this case, is finite (con^esponding to its value at the pole) and larger than r. In all 
panels we had set r = 10~^s~^; the results of ten realizations of the numerical solution (empty symbols) 
are shown for each value of T or fcwind- The vertical dashed line in panels (a), (b) indicates the time period 
associated with the Coriolis frequency at the pole. 
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FIG. 3: (a) Comparison between the second moment of the current amphtude as calculated analytically 
(solid line), numerically (red circles) and by the MITgcm model (blue triangles) in a simple lake. The wind 
stress amplitude was drawn from a Weibull distribution with k = 1.5. (b) The Weibull k parameter of the 
current amplitude distribution versus the k parameter of the wind stress for different coiTclation times. In 
both panels the Coriolis parameter, /, corresponds to its value at lat = 30° and r = 10~^s~^. 



linearity of the model, the current components have Gaussian distributions which are fully char- 
acterized by the mean and variance. Note that the result of eq. (fTTI) holds for any distribution of 
the wind stress components as long as the two-point correlation function is given by eq. It is 
clear, similar to the first idealized case discussed above, that the second moment vanishes at the 
limit of a short wind stress correlation time, (7 ^ r and 7 ^ /). 

Test on detailed GCM model: We have performed two types of numerical simulations to vali- 
date the analytical derivations presented above. The first type was a simple integration of eqs. 
As for the second type, we have used a state-of-the-art oceanic GCM, the MITgcm H, to test 
the applicability of our analytical derivations when the spatial variability of bottom topography 
and the nonlinearity of the ocean dynamics 111 ill are taken into account. We have considered a 2D 
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parabolic basin with a maximum depth of 90m that is situated in a plateau of a depth of 100m. 
The boundaries of the domain are open, and a grid with 50 x 50 points and a resolution of 1km is 
assigned. The integration time step is 10s, and the overall integration time is two years. As shown 
in the figures, these simulations exhibit close similarity with the analytical derivations, suggesting 
that these may also be relevant to more realistic cases, in the presence of bottom topography and 
nonlinearities. We have also performed similar simulations with closed boundaries and with a 10 
times coarser spatial resolution; we obtained similar results in the center of the domain but not 
close to the periphery, suggesting that our derivations are valid in the open ocean but not close to 
shores. 

Discussion and Summary: The results presented above highlight the important role that the tem- 
poral correlations of the wind stress has on the statistics of surface ocean currents. The behavior 
at the limits of long-range temporal correlations (rT ^1) and short-range temporal correlations 
(T —J- 0) are intuitive, once derived. The existence of an optimal correlation time, at which the 
average current amplitude is maximal, is less trivial. 

In Fig. [2ta), we show the second moment of the current amplitude for the isotropic, step-like wind 
stress model. The analytical results (eq. ([7])) are compared with the numerical solution of the 
Ekman model (eq. dH)) and MITgcm modeling of the currents in a simple artificial lake (details of 
which are provided above). One can see that the dependence of the average current amplitude on 
the constant wind stress duration is nonmonotonic (due to the fact that the Coriolis effect is sig- 
nificant - I/I > r) and that there is an excellent agreement between all the results. At the equator, 
where / = 0, the second moment increases monotonically to (Irl^) / {r^pD as a function of T. 
It was previously argued that, under certain conditions, the wind amplitude (directly related to 
the wind stress) PDF is well approximated by the WeibuU distribution [3]; we thus choose, in 
our demonstrations, a WeibuU distribution of wind stress. In Fig. [2tb), we present the WeibuU k 
parameter ll^ of the current distribution versus the constant wind stress duration for two different 
values (k = 1 and A; = 2) of the wind stress WeibuU k parameter (again, we only present here the 
results of the isotropic case). One can see that, for short correlation times, the current amplitude 
exhibits a Rayleigh distribution (k = 2), independent of the wind stress distribution. This corrre- 
sponds to Gaussian distributions of the current components. On the other limit of long constant 
wind stress periods, the current distribution WeibuU k parameter converges to that of the wind 
stress (not shown for A; = 2). 

The above results were obtained for the maximal value of the Coriolis parameter, i.e., / at the pole. 
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In Fig. [2tc,d) we show the WeibuU k parameter of the current amplitude distribution versus the 
WeibuU k parameter of the wind stress for different values of the constant wind stress duration, T, 
and the two limiting values of the Coriolis parameter, /, at the equator and at the pole. Here again, 
one observes an excellent agreement between the numerical and the analytical results for both 
values of /. The limits of short and long temporal correlations of the wind at which /^current = 2 
and /^current = ^wind are clearly demonstrated. 

In Fig. [3l we present similar results for the case of an exponentially decaying correlation function 
of the wind stress. In Fig. Oa), we present the second moment of the current amplitude versus the 
decay rate of the correlation function. The existence of an optimal decay rate, for which the aver- 
age current amplitude is maximal, is demonstrated for this case as well. One notable difference is 
the absence of the wiggles which appeared in the step-like model. For this case, one may easily 
find that the maximal average amplitude of the currents is obtained for 7 = |/| — r, assuming that 
I /I > r. In Fig. Ob), we present the WeibuU k parameter of the current amplitude distribution 
versus that of the wind stress distribution. Here again, we obtain an excellent agreement between 
the predicted and numerically obtained limiting behaviors. 

In summary, we have shown, in two idealized cases, that the PDF of wind-driven ocean currents 
depends on the temporal correlations of the wind. For short-range correlations, the current speed 
approaches zero, and the PDF of its components is Gaussian. For long-range temporal correlations 
of the wind, the currents' PDF is proportional to the wind-stress PDF. The model we used highly 
simplifies real ocean dynamics, yet our numerical results suggest that, qualitatively, the above 
conclusions should also be valid in more general, realistic cases. Analysis that is based on the 
depth dependent model is a natural extension of the current study and will allow us to compare the 
analytical results to altimetry based surface currents. 
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